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Investigating different structures of the ^^(lOGlO) and Zi,{10650) 

Chun-Yu Cui, Yong-Lu Liu and Ming-Qiu Huang 

Department of Physics, National University of Defense Technology, Hunan 410073, China 

The recently observed narrow resonance Zb( 10610) is examined with the assump- 
tions both as a B*B molecular state and a [6d][6?2] tetraquark state with quantum 
numbers I^ J^ = 1+1+. Possible interpolating currents are constructed to describe 
the Zfe(10650) as an axial-vector B*B* molecular state or a [fed] [fen] tetraquark state. 



^D ■ Using QCD sum rules (QCDSR), we consider contributions up to dimension six in the 

operator product expansion (OPE) at the leading order in ag- The mass is obtained 



as (10.44±0.23) GeV for molecular state and (10.50±0.19) GeV for tetraquark state, 
both of which coincide with the ^^(10610). The results m^,g, = (10.45 it 0.31) GeV 
and wirj^irj^i = (10.48 it 0.33) GeV are consistent with the ^^(10650). 

PACS numbers: 11.55.Hx, 12.38.Lg, 12.39.Mk 



^ . I. INTRODUCTION 



^ I In the past years, the Babar, Belle, CLEO, DO, CDF and FOCUS collaborations re- 

cn . ported many charmonium-like states which stimulate theorists' interests in revealing their 

underlying structures [l|. There comes to the consensus that these new states are not 
regular mesons or baryons. Various approaches have implied that a possible candidate is 
the exotic state, which means a complex structure such as a molecule, a tetraquark or 
a hybrid. This statement indicates that the counterpart may exist in the bottomonium 
region. 
1^ After the two charged bottomonium-like resonances Zft(10610) and ^^(10650) being ob- 

;_, I served in the 7r^T(nS') (n = 1,2,3) and n^hi^irnP) [m = 1,2) mass spectra by Belle 



Collaboration |2|, many atternpts have been made to investigate their possible configura- 



tions with various models |3l-ll3l|. most of which support the B*B^*^ molecular structure 
with J^ = 1+. In ref. [3|, an explanation of two charged bottomonium-like resonances 
in terms of cusps at B*B^*^ channel is presented. In Ref. [4], the authors discuss the 
special decay behaviour of the J=l S-wave B*B^*^ molecular states and study radiative 
transitions from T(55') to molecular bottomonium based on the heavy quark symmetry. 
In Ref. |6|, Zb(10610) and ^^(10650) are interpreted as tetraquark states in the framework 
of chromomagnetic interaction Hamiltonian model. In Ref. jl2|, the authors investigate 
the mass spectra of S-wave bound states BB* and B*B* systems with quantum numbers 
I{J^'") = 1(1^~) in the framework of a chiral quark model. All in all, quantum numbers 
compatible with the experiment are the fundamental ingredients in theoretical analysis of 
composite particles within any specific model. We notice that Zf,(10610) was observed in 



the T(55') —7- ^^,(10610) + vr" decay process. Since the T{5S) has negative G-parity, due 
to emission of the pion, the newly observed ^^(10610) has positive G-parity. In consider- 
ation that it is charged, the known quantum numbers of ^^(lOGlO) are J*^ = 1+. Angular 
distributions analysis favors the J^ = 1+ assignment for Zf,(10610) |[2|]. The ^6(10610) 
considered as a B*B molecular state in Ref. |7| has J^ = l"*" with no definite quantum 
numbers of isospin and G- conjugation. Considering all their known quantum numbers 
I^J^ = 1^1^, we construct a molecular state with a B and a B* mesons using the config- 
uration B^B*^ + B*~B^. It is difficult to construct a suitable axial-vector style molecular 
current interpolating the state Zfe(10650) using both B* and B* fields. A possible interpo- 
lator is supposed to describe the axial-vector style molecular states B*B*: 

where qi stands for light quarks. Performing the parity transformation to the current, it 
satisfies the condition Pj^[x)P~^ = j^{x). 

Mass property is expected to be helpful for understanding the configuration of 
bottomonium-like resonances Z^. In the hadronic scale, it is difficult to get reliable the- 
oretical estimate for the mass using the perturbative QCD. Therefore, we need some 
non-perturative methods to describe the non-perturative phenomena. QCDSR 
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18| is powerful since it is based on the fundamental QCD Lagrangian. We notice 



that in the case of X(3872) with quantum numbers J^*" = 1"J^ the authors study 



its mass with QCDSR using the configuration D^D*^ — D*^D^ 19| and the configura- 
tion [eg] 5=0 [eg] 5=1 -|- [eg] 5=1 [eg] 5=0 |20|. Different from above situation, we consider the 
Zf,(10610) with known quantum numbers I^J^ = 1+1+ with no definite charge conjuga- 
tion. This work is devoted to investigate the masses of ^^(10610) and Zf,(10650) in the 
QCDSR, both a BB* molecular state and a [6'u][M] tetraquark state being assumed. 

The rest of the paper is organized as three parts. The QCDSR for the Zb{10610) is 
derived in Sec. [TTl with contributions up to dimension six in the OPE. The numerical 
analysis is presented to extract the hadronic mass at the end of this section. Sec. IIIII is 
organized for the QCD sum rules of the Zf,(10650), with numerical estimation of the mass. 
Sec. IIVI is the summary and conclusion. 

II. QCD SUM RULES FOR ^^(lOeiO) 
A. molecular state QCD sum rules 



In the previous work |7|, Zb(10610) has been studied as a B*B molecule with J^ = 1"*" 
in the framework of QCDSR. However, the information of isospin and G-parity, which have 
been confirmed by the experiment, is omitted in constructing the interpolating current. In 
this subsection, the Zf,(10610) resonance is considered by constructing a proper interpolator 



as a B*B molecule with all the known quantum numbers I J = 1"*"1+ [2|]: 

f = i=[(m756)(67^d) + (M,7A)(&fcn'4)]. (2) 

In the QCDSR approach, the mass of the particle can be determined by considering the 
two-point correlation function 

U^'^iq') =zJd'xe''^-^0\T[f{x)f+iOm- (3) 

Lorentz covariance implies that the two-point correlation function can be generally param- 
eterized as 

n^^(g^) = (^ - 9n^^'\q') + ^nW(g2). (4) 

qZ gZ 

The term proportional to Qf^,^ will be chosen to extract the mass sum rule, since it gets 
contributions only from the l"*" state. The QCD sum rule attempts to link the hadron 
phenomenology with the interactions of quarks and gluons. It contains three main ingredi- 
ents: an approximate description of the correlation function in terms of intermediate states 
through the dispersion relation, a description of the same correlation function in terms of 
QCD degrees of freedom via an OPE, and a procedure for matching these two descriptions 
and extracting the parameters that characterize the hadronic state of interest. 

In the phenomenological side, the correlation function is calculated by inserting a com- 
plete set of intermediate states. Parameterizing the coupling of the Z state to the current 



j^ as 





{Oj,Z) = Xe,. 


•^^^(5^) can be expressed as 




n«(g^) = 


A2 I r^ Imn(i)Phen(s) 


Ml-q^ ' rcJso s-q^ 



(5) 



(6) 

where Mz denotes the mass of the molecular state, and sq is the threshold parameter. 
In the OPE side, 'US^\q^) can be written as 

where the spectral density is p'^^^{s) = -Imli^ ){s). Making quark-hadron duality as- 
sumption and a Borel transformation, the sum rule is obtained by matching the two sides: 

JAml 

with M^ the Borel parameter. 



In the OPE side, we work at the leading order in a^ and consider vacuum condensates 
up to dimension six, with the similar techniques in Refs. [2l|. In order to consider the 
isospin violation, we keep the terms which are linear in the light-quark masses rriu and 
171^. After some tedious OPE calculations, the concrete forms of spectral densities can be 
derived: 

with 

Z TT ^ OLmin Cx 'J Pmin P 



+ TTm^ I T7 L ^5j(a + /3- l)(m„a + md/3 
+ mua(3 + mda(3 + Sniua + 3md(3)r{mb, s)^ 



(^^)f.^ - JML^. r-^ /'""^r.,^mn^.,^«w^, .^^ 



p'''Ks) = -l^^b / -^ / -^(a + /3)(1 + a + P)rim,, s) 



+ l^^rr("^« + ^d) / -7- -ri^Q - a(l - a)s 

r-K^ Ja^i^ a\v — a) 

3(qq) , , r^rna:c da f^-°' dl3 , ,2 

'niu + md) / — / — -r(mfe,s) , 



Z 7T Jamin a J Pmin P 



Z TT Joimin Oi 

3{gqcT-Gq) r^ma. fi-^ dl3 
+ ^ITi ^b da — r(mb,s), 

"^ ^^n "^M ^" / -j^{(^ + P)r{mb,s} 



Z TT '' Olmin •'Pmin, P 

(a^G^) f^rnax fl-Oi dB 

+ ini ^«/ -|(l-«-/3)(l + « + /5)r(m,,s), (10) 

Z TT '' Ctmin •' Pm.in P 

with r{mh,s) = {a + /3)m^ — a/3s. The integration limits are given by ttmm = (l ~ 

y'l -4m^/s)/2, a^ax = (l + y'l -4mf/s)/2, and /3mj„ = aml/{sa - ml). 

Taking the derivative of Eq.( 1T7|) with respect to -^ and then dividing by itself, we arrive 



TABLE I: Upper limits in the Borel window for the I^J^ = 1+1+ B*B* current obtained from 
the sum rule for different values of ^/sQ. 



V^ (GeV) 


MLAGeM') 


11.0 


8.9 


11.1 


9.2 


11.2 


9.6 


11.3 


9.9 


11.4 


10.3 


11.5 


10.7 



at the mass of the molecular state 



Ml 



so 



4m? 

b 



dsp' 



OPE 



se 



-s/M^ 



so 



b 



dsp 



OPE -s/M^ 



:iii 



Before the numerical analysis of the equation f lTT|) . we first specify the input parameters. 
The quark masses are taken as rriu = 2.3 MeV, rrid = 6.4 MeV, and rrih = (4.24 ±0.06) GeV 
22|. The condensates are (uu) = (dd) = (qq) = -(0.23±0.03)3 GeVf , (gqcr-Gq) = ml (qq), 
ml = 0.8 GeV^ {g^G^) = 0.88 GeV^ and {g^G^) = 0.045 GeV^ [16]. Complying with the 
standard procedure of the sum rule, the threshold sq and Borel parameter M^ are varied to 
find the optimal stability window. There are two criteria (pole dominance and convergence 
of the OPE) for choosing the Borel parameter M^ and threshold Sq. 

The contributions from the high dimension vacuum condensates in the OPE are shown 
in FigJH We have used sq > 121 GeV^. From this figure it can be seen that for M^ > 
8.10 GeV^, the contribution of the dimension-6 condensate is less than 5% of the total 
contribution and the contribution of the dimension-5 condensate is less than 20% of the 
total contribution, which indicate a good Borel convergence. Therefore, we fix the uniform 
lower value of M^ in the sum rule window as M^j„ = 8.10 GeV^. The upper limit of M^ 
is determined by imposing that the pole contribution should be larger than continuum 
contribution. Figj2] shows that the contributions from the pole terms with variation of the 
Borel parameter M^. We show in Table |T]the values of M^^^. for several values of ^/so■ 
In Figjni we show the molecular state mass, for different values of ^/so, in the relevant 
sum rule window. It can be seen that the mass is stable in the Borel window with the 
corresponding threshold a/sq. The final estimate of the I^ J^ = 1+1+ molecular state is 
obtained as 



Mz = (10.44 ± 0.23) GeV. 



(12) 
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FIG. 1: The OPE convergence for the I J = l"*"!^ molecular state. The contributions from 
different terms with variation of the Borel parameter M^ in the OPE. The A and B correspond to 
the contributions from the D = 5 term and the D = 6 term, respectively. The notations a, /3, 7, 
A, p and a correspond to the threshold parameters sq = 121.00 GeV^ 123.21 GeV^ 125.44 GeV^, 
127.69 GeV^, 129.96 GeV^ and 132.25 GeV^, respectively. 
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FIG. 2: The contributions from the pole terms with variation of the Borel parameter M^ in the 
case of molecular state. The notations q, /3, 7, A, p and a correspond to the threshold param- 
eters So = 121.00 GeV^ 123.21 GeV^, 125.44 GeV^, 127.69 GeV^, 129.96 GeV^ and 132.25 GeV^ 
respectively. 



B. tetraquark state QCD sum rules 



As pointed out in ref. |6i], if the charmonium-like tetraquark X(3872) and Z+(4430) are 
really tetraquark states, the bottomonium-like tetraquark should exist. A possible current 




M^IGeV^] 



FIG. 3: The mass of the I J = V^V^ molecular state as a function of M^ from sum rule (jlip . 
The notations a, /3, 7, A, p and cr correspond to the threshold parameters sq = 121.00 GeV , 
123.21 GeV^ 125.44 GeV^ 127.69 GeV^ 129.96 GeV^ and 132.25 GeV^ respectively. 



interpolating a I^ J^ = 1+1+ tetraquark state with [6(i]5=o and [6m]5=i fields is given by 
J, = '-^[{d^aC^,h){u,^,Cb^) - {d:C^,h){u,^,Cb^)] , (13) 

The spectral densities of the sum rule are obtained with the same standard procedure of 
the approach as what was done in subsection III Al for the molecular current: 



P^^^(^) 



P 



pert. 



s) + p^^'^^s) + p 



(g^G^), 



s)+P 



(gqa-Gq) , 



s)+p<«"«>'(s)+p 



{9'G^), 



(14) 



with 



pPe^H^) 



"^ da /■!-" d/3 



{l~a-(3){l + a + (3)r{mb,sY 



9107^6 r n/S a R3 

nib /■"""'" da /■!-" d(3 , n .x, 1 ni 

Z TT J (Xmin a J Pmin P 

+ m„a/3 + mdal3 + Sm^^a + 3mdl3)r{mb, s)^, 



p^^^^s) 






da 



afl — a 



-[ml-a{l-a)sf 



(qq) r°'^o,^ da f'^-°' df3 , ^,,^ ^, , ,n 

- ' -T -^(a + /3)(l + a + /3)r(m6,s)2 






p<^'^'>(.) 



267r4 

{g'G') 



min 



da /"i-" dp , , 

— / -n-r[mb,s) 

\jL ^ Pmin f 



mu + rrid) I ^^ / ^r(mfc, s)^ 



m. 



3*297r6 ^ 



^min ^ Pmin /-^ 

/•I— a (^fl 

(ia / — (1 - a - /3)(1 + a + /3)r(mb, s) 



+ 



3 * 21V 



12^6 '"b 



amax da p-"' dl3 , ^., ^, , , 

m,, / — / -^[1 - a - p)[3 + a + p)r[mb,s) 

in LX J p^i 



P' 



TABLE II: Upper limits in the Borel window for the I^J^ = 1+1^" tetraquark current obtained 
from the sum rule for different values of ^/s(). 



V^(GeV) 


MLAGeY') 


11.0 


8.7 


11.1 


9.1 


11.2 


9.4 


11.3 


9.8 


11.4 


10.2 


11.5 


10.6 



p(9g--Gg>^3) 



p<'^^>'(s) 



p^^'^^'^s) 



+ 



3 * 2i07r6 Ja 
{gqa ■ G q) 

2%4 

{gqa-Gq) 
+ ^TT^ rUb 



<^^ da /"i-" d(5 



^ ^ Pmin f~^ 



{2a + 2P-l)r{mb,sy 



-nib 



Ol-max dot 9 

— [m^ — a(l — a)s\ 

-m.in. '-^ 



2H^ 



amax da r^-°' 



{gqa ■ Gq) 



267r4 



rrib 



min ^ ''Pmin 



da 



mm 



Pmin /-^ 



df3r{mb, s) 

a + l3)r{mb,s) 



{gqa ■ Gq ) 
3 * 2%^ 
{qq)^ 



-nib 



amax (l(y rl—a 

— dP{a + P + l)r{nib,s), 

min ^ "^ Pmin 



3 * 227r2 

{g'G') 

3 * 2i07r6 

{g'G') 



m^Jl-Aml/s, 



m. 



+ 



da 



da 



l-a 



Pr, 



^«(l-«-/3)(l + « + /?) 



3*2 TX Jamin J Pmin P 



1-° dB 

-^(l-a-/3)(l + a + /3)r(mb,s). (15) 



The mass sum rule is obtained as the same form of Eq. fITT]) . In this case, from FigJUwe 
see that we obtain a reasonable OPE convergence for M^^^ = 8.90 GeV^. From this figure 
it can be seen that for M^ > 8.90 GeV^, the contribution of the dimension-6 condensate is 
less than 9% of the total contribution and the contribution of the dimension-5 condensate 
is less than 25% of the total contribution. FigjS] shows that the contributions from the 
pole terms with variation of the Borel parameter M^. The upper limits of M^ for each 
value of ^/sq are given in Table HTl The table also indicates that ^/s^ > 11.1 GeV to ensure 
^max — ^min- ^^ FigJHl wc show the mass Mz depending on the Borel mass for several 
threshold values ^/sQ. It can be seen that we get a good Borel stability for Mz- The 
numerical result is 



Mz = (10.5 ±0.19) GeV. 



(16) 
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FIG. 4: The OPE convergence for the I J = l"*"!"*" tetraquark state. The contributions from 
different terms with variation of the Borel parameter M^ in the OPE. The A and B correspond to 
the contributions from the D = 5 term and the D = 6 term, respectively. The notations a, /3, 7, 
A, p and a correspond to the threshold parameters sq = 121.00 GeV^, 123.21 GeV^, 125.44 GeV^, 
127.69 GeV^, 129.96 GeV^ and 132.25 GeV^, respectively. 
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FIG. 5: The contributions from the pole terms with variation of the Borel parameter M^ in the 
case of tetraquark state. The notations a, /3, 7, A, p and a correspond to the threshold param- 
eters So = 121.00 GeV^ 123.21 GeV^, 125.44 GeV^, 127.69 GeV^, 129.96 GeV^ and 132.25 GeV^ 
respectively. 



It is noticed that the result is very close to the one in Eq. (fT2|) . The reasons lies in 
the following arguments: The mass sum rule is expressed as a division of two parts as 
Eq.fITT]). However, in the two interpolating currents involved in the paper, the dominant 
contributions to the sum rules are the perturbative term and D = 3 condensate term, 
which are proportional by 3/4 that will be canceled in the process of division for both 
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M^IGeV^]" 



FIG. 6: The mass of the tetraquark state as a function of M^. The notations a, /3, 7, A, 
p and a correspond to the threshold parameters sq = 121.00 GeV^, 123.21 GeV^, 125.44 GeV^, 
127.69 GeV^ 129.96 GeV^ and 132.25 GeV^ respectively. 

cases. The deviation comes from the D = 4 and D = 5 condensates terms which play 
subdominant roles in the result. It tacitly suggests that two-point sum rules is unable to 
give explicit distinction whether the exotic is a molecular state or a tetraquark state. 



C. General analysis on the mass sum rule for a four-quark state 

From the above study, we notice that for a four-quark state, there are two main candi- 
dates for its possible configuration which are a molecular state and a tetraquark state. In 
consideration of the requirement of the quantum numbers and the color singlet, the general 
form can be written as qiTiQQTjq2 for a molecular state and eabc^dec{(li^CTiQ^){q2TjCQ'^^) 
for a tetraquark. In the expressions ^i{j) is a general form of a gamma matrix, and eabc(def) 
is the usual antisymmetric three order tensor with a, 6, c, rf, e being color indices. 

For a concrete form, there are ten possible combinations of '^i(j) for quantum numbers 
J^ = 0~, J^ = 0"*", J^ = 1~ and J^ = 1"*". The ten possible currents for the two 
configurations are shown in Tab JIIIi To get the sum rule, both the dispersion relationship 
and the quark-hadron duality approximation are used. The phenomenological part can be 
expressed as the integral of the spectral density: 



;^2g-A/|/M2 



so 



Ami 



dsp 



OPE, 



s]e 



-s/M'^ 



(17) 



After some tedious OPE calculations, the concrete forms of spectral densities can be 
derived: 



pert 



p^^^(s)=pP«^(s)+p^««>(s)+p' 



(ot>( 



higherorder , 



S). 



:i8) 
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TABLE III: Possible combinations of ^iij) for different J^ quantum numbers. In each row, 
qiTiQQTjq2 stands for a molecular state and eabc^dec{Qi^CTiQ''){q2^jCQ^'^) stands for a 
tetraquark state. 



J^ 


li 


7i 


JP = 0- 


75 


1 


jp = o- 


7/. 


7/^75 


JP = 0+ 


75 


75 


JP = 0+ 


1 


1 


JP = 0+ 


7/. 


7m 


JP = 0+ 


7^75 


7/^75 


JP = l- 


1 


7m 


JP = 1- 


75 


7m 75 


JP = 1+ 


75 


7m 


JP = 1+ 


1 


7m75 



To get the mass sum rule, we need to take the derivative of Eq.([T 
and then divide by itself. The final result is expressed as follows: 



Ml 



SQ 



Avaj- 



dsp' 



OPE 



se 



-s/R4'' 



so 



dsp 



OPE -s/M'^ 



with respect to -gj 



(19) 



A simple calculation shows that the two dominant contributions p^ {s) and p^^'^^s) 
are same for the two configurations except for a factor. As in the analysis of the sum rules, 
Borel transformation is used to suppress both the higher resonance contribution and make 
the OPE convergence better. One of the criteria for the choice of the Borel parameter is 
that in the working region the higher order contributions of the vacuum condensates are 
required to be less than about 20%. However, the errors of the sum rule come from the 
variation of the Borel parameter, the threshold and the input parameter, which may bring 
an uncertainty bigger than 20%. Thus, deviation from second order vacuum condensate 
contributions can not give a definite judge on the fact which configurations can be in 
accordance with the exotic state concerned. 

The two different configurations are quite different and are governed by different dy- 
namics. Unfortunately, the mass sum rules cannot give us the useful detailed information 
on this deviation. Therefore, thorough understanding of the internal structure of a exotic 
state in the framework of QCDSR requires to investigate various decay processes that may 
contain more detailed information of the state. 
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III. QCD SUM RULES FOR Zfc(10650) 

A. molecular state QCD sum rules 

Due to the mass of ^^(10650) lies close to the B*B* threshold, B*B* molecule seems 
to be a natural candidate. Based on B*B* molecular type structure, ZftflOGSO) is studied 
within one-boson-exchange model and XEFT theory in Ref. [4, llO, Il3|. In order to inves- 



tigate the state of this configuration with QCDSR, we construct a possible interpolator: 

f[x) = e^'"'^{u{x)t^Mx))D^{b{xh^d{x)). (20) 

Following the standard procedure of QCDSR, we obtain the concrete forms of spectral 
densities as follows: 

pO^^(s) = pPert (,) + p(5a>(,) + p(^^G^>(,) + p(5?-^'?)(s) + p(^'^)\s) + p<^'^'>(s), (21) 
with 



5*2 Vr Jamin O! J I3min P 



2°7r^ "'a™,,, a 1 — ar 



+ ^mj — / ^(«/3 + /32-2a-2/3)r(mfe,s)^ 

O * Z 71 Jamin Oi J (5,nin P 

- 2a^ + 7af3^ + 4/3^ + 8(3^)r{mb, s) 

{a'^G'^) 1 r""''' da f^-°' d(3 , f, ^ .^ ^ . .0 4^ 
^ -"^6 / — r / -^{a^ + 3a'^P -6a^ -18a^f3 



0*2 vr Jamin Oi J Pmin P 

+ 21a^ + lOa^/3^ - Qa^/J^ + 15a^/3 - lOa^ + I2a^[3^)r{mj,, sf 
+ Tr^/ — / -g(2«' + 17«'-4a/3)r(m,,s)3, 

+ 2%4 "^^L„ ^^-^a.K - «(1 - «).] 



097^4 ° / ^2 /o /32 

3 * 257r2 



P^^"'^^'(^) = -^^^,{^< + rnls)^^-^</s, 



P^^'^'H-) = -^^^,mir^^dar"^{7a' + 27a'p-39a' + 21a(3' 
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TABLE IV: Upper limits in the Borel window for the axial- vector B*B* molecular current ob- 
tained from the sum rule for different values of v^- 



B*B* 

Vi^(GeV) 


B*B* 


11.0 


8.5 


11.1 


8.8 


11.2 


9.1 


11.3 


9.4 


11.4 


9.7 



- 30a/3 + 9a + f3^- 3/3^ + 3(3- l)r{mb, s) 



+ 






da 
a 



l-a 



dp_ 



(l-a-/3)V(mb,s)' 



m\ 



The contributions from the high dimension vacuum condensates in the OPE are shown 
in FigJTl We have used sq > 121.00 GeV^ in the B*B* channel. From this figure it can be 
seen that for M^ > 8.0 GeV^ in the B*B* channel, the contribution of the dimension-6 con- 
densate is less than 7% of the total contribution and the contribution of the dimension-5 
condensate is less than 18% of the total contribution, which indicate a good OPE con- 
vergence. Therefore, we fix the uniform lower value of M^ in the sum rule window as 
^rain — 8.0 GeV^ in the B*B* channel. FigJH] shows that the contributions from the pole 
terms with variation of the Borel parameter M^ and the threshold parameters sq. The 
upper limit M^^^ is determined with the requirement that the pole term contribution is 
bigger than 50% of the total contribution. We show in Table [IV] the values of M^^^ for 
several values of ^/sq in the B*B*. In FigJHl we show the molecular state masses, for differ- 
ent values of ^/so, in the relevant sum rule window. It can be seen that the mass is stable 
in the Borel window with the corresponding threshold ^/so■ 

Up to now we have fixed the values of the quark masses and condensates. To make the 
results more reliable, we also consider errors from the uncertainties of input parameters mi,, 
(qq), and ml. Taking into account both uncertainties of input parameters and uncertainties 
due to the continuum threshold parameter and Borel window, we finally arrive at 



M, 



B*B* 



(10.45 ±0.31) GeV. 



(23) 



B. tetraquark state QCD sum rules 



As the case of Zf,(10610), there are also two possible configurations for ^^(10650). In 
this subsection, we consider the tetraquark structure, which can be described by the inter- 
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M2[Gev2] 



(a) 



(b) 



FIG. 7: The OPE convergence for the axial-vector B*B* molecular states. The contributions from 
different terms with variation of the Borel parameter M^ in the OPE. The A and B correspond 
to the contributions from the D = 5 term and the D = 6 term, respectively. The notations a, /3, 
7, A and p correspond to the threshold parameters sq = 121.00 GeV^ 123.21 GeV^ 125.44 GeV^ 
127.69 GeV^ and 129.96 GeV^, respectively. 
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FIG. 8: The contributions from the pole terms with variation of the Borel parameter M^. The 
notations a, /?, 7, A and p correspond to the threshold parameters sq = 121.00 GeV^, 123.21 GeV^, 
125.44 GeV^ 127.69 GeV^ and 129.96 GeV^, respectively. 



polator: 



fix) 



e^'"'^'-^^^^{u{x)lC^X^),)DM^)dlpCh{x) 



V2 



(24) 



After some tedious calculations, the concrete forms of spectral densities can be derived: 



^OPE( 



.pert. 



s) + p'^^''\s) + p 



{g^G'^), 



s)+P 



{gqa-Gq) , 



s) + p^^'i^\s) + p 



{g'G^), 



(25) 



with 



O * Z 71 JOmin W J Pmin P 
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FIG. 9: The masses of the B*B* molecular state as a function of M^. The notations a, /3, 7, 
A and p correspond to the threshold parameters sq = 121.00 GeV^ 123.21 GeV^, 125.44 GeV^ 
127.69 GeV^ and 129.96 GeV^. 
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-(l-a-/3)V(mfe,s)' 
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(26) 

It can be found that in the two configurations involved for ^{,(10650), the dominant 
contributions to the sum rules are the perturbative term and D = ?> condensate term, which 
are proportional by 3/4 that will be canceled in the process of division for both cases. The 
deviation comes from the D = A and D = 5 condensate terms which play sub dominant 
roles in the result. It tacitly suggests that the same conclusion can be drawn as in SecHTl 
that two-point sum rules is unable to distinguish whether Zfe( 10650) is a molecular state 
or a tetraquark state. Our final numerical result is 

M[,,][-,^] = (10.48 ± 0.33) GeV. (27) 

IV. SUMMARY AND CONCLUSION 

By assuming Zf,(10610) as both a B*B molecular state and a [6(i][6u] tetraquark state 
with quantum numbers I'^J^ = 1^1^, the QCDSR approach has been applied to calculate 
the mass of the resonance. Our numerical results are Mz = (10.44 ± 0.23) GeV for 
molecular state and Mz = (10.50 ± 0.19) GeV for tetraquark state. Both of the results 
are compatible with the experimental data of ^(,(10610) by Belle Collaboration. We also 
construct possible interpolators to describe the Zf,(10650) as both an axial-vector B*B* 
molecular state and an axial- vector [6(i][6n] tetraquark state. Our numerical result are 
Mb*b, = (10.45 ±0.31) GeV and m[f,d][bu] = (10.48 ±0.33) GeV, which are compatible with 
the experimental data of Zft( 10650). 

The calculations indicate that the mass sum rule could not distinguish Zf,(10610) 
(^5(10650)) between a B*B {B*B*) molecular state and a [M][6'u] tetraquark state. The 
clarification of the configuration of these two states requires further analysis on the decay 
channels Z6(10610)(Zfc(10650)) -^ 7i^T{nS) (n = 1,2,3) and Zfe(10610)(Z6(10650)) -^ 
TT^hb{mP) {m = 1,2), which contain more detailed dynamical information. 
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